Noise Enhanced Stability in Fluctuating Metastable States 
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We derive general equations for the nonlinear relaxation time of Brownian diffusion in randomly 
switching potential with a sink. For piece-wise linear dichotomously fluctuating potential with 
metastable state, we obtain the exact average lifetime as a function of the potential parameters and 
the noise intensity. Our result is valid for arbitrary white noise intensity and for arbitrary fluctuation 
rate of the potential. We find noise enhanced stability phenomenon in the system investigated: the 
average lifetime of the metastable state is greater than the time obtained in the absence of additive 
white noise. We obtain the parameter region of the fluctuating potential where the effect can be 
observed. The system investigated also exhibits a maximum of the lifetime as a function of the 
fluctuation rate of the potential. 

PACS numbers: 05.40.-a,02.50.-r,05.10.Gg 
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I. INTRODUCTION 



Activated-escape in systems with metastable states un- 
derlies many physical, chemical and biological problems. 
Examples are crystal growth, tunnel diode, lasers, quan- 
tum liquids, spin systems, protein folding and polymer 
physics P, Q ■ The most interesting and stubborn case of 
metastable state is one described by time dependent po- 
tential, which fluctuates on a characteristic time scale 
that may vary over a large range. In particular, the 
metastable states with fluctuatingbarriers are common 
to chemical and biological models 0, Q , and to a wide 
range of physical problems, such as nonequilibrium trans- 
port models, molecular dissociation in strongly coupled 
chemical systems Jal , ratchet models for the action of 
molecular motors [3 , noise in microstructures and gener- 
ation process of ca"ier traps in semiconductors Siia. 
The escape from metastable state with fluctuating or ran- 
domly switching barrier was studied in the past mainly 
by well known mean first passage time (MFPT) tech- 
nique. In Refs. [H 113 exact results for the MFPT of 
escape process over fluctuating barrier that switches be- 
tween two configurations have been obtained. However, 
the MFPT method requires the implication of absorbing 
boundary in the system, and it does not take into account 
the inverse probability current through this boundary. 
The nonlinear relaxation time (NLRT) method is devoid 
of this disadvantage ■ Nevertheless the theory for the 
NLRT is not well developed and the equations for the 
NLRT are unknown for the case of time varying poten- 
tial. 

In the present paper we derive general equations for 
the NLRT for potentials randomly switching between two 
arbitrary configurations with a sink. Wc find the exact 
solution of these equations for a piece-wise linear poten- 
tial flipping between unstable and metastable configu- 
rations, for arbitrary white noise intensity and fluctua- 
tion rate of the potential. Analyzing this exact result 



we focus on the noise enhanced stability (NES) effect, 
which implies that the system remains in the metastable 
state for a longer time than in the absence of additive 
white noise, and the lifetime of the metastable state 
has a maximum at some noise intensity. This effect, 
which cannot be described by Kramers-like behavior, was 
observed and investigated theoretically and experimen- 
tally in various physical systems and mainly concerning 
MFPT in periodically or randomly driven metastable 

states 0, la El III 111 13 13 111 ES IM m, 23, 2|. 

In these papers the nonmonotonic behavior of the av- 
erage escape time was observed: (i) in physical systems, 
like tunnel diode 0|, and Josephson junction ^3, where 
the influence of thermal fluctuations on the superconduc- 
tive state lifetime and the turn-on delay time for a single 
Josephson element with high damping was investigated; 
(ii) in chemical systems, like the one-dimensional return 
map of the Belousov-Zhabotinsky reaction, by investi- 
gating the behavior of the lenght of the laminar region 
as a function of the noise intensity and (iii) in bio- 
logically motivated models, such that investigated in ref 
|3j, where the overdamped motion of a Brownian parti- 
cle moving in an asymmetric fluctuating potential shows 
noise induced stability. 

Here we study the NES phenomenon for the NLRT 
in randomly switching metastable state, and we obtain 
analytically the region of system parameters, where this 
effect takes place. We find also resonant activation phe- 
nomenon by investigating the mean lifetime as a func- 
tion of switchings mean rate. Moreover, we find that 
the NLRT exhibits a maximum as a function of barrier 
switching rate. This new resonant-like phenomenon is 
related to the NES effect |l|l3. 

The paper is organised as follows. In the second section 
we derive the general equations for the nonlinear relax- 
ation time of Brownian diffusion in randomly switching 
potential with a metastable state. In the third section we 
analytically derive the mean lifetime for piece- wise linear 
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potential. In the next section we obtain the condition to 
observe the NES phenomenon and investigate the behav- 
ior of the mean hfe time as a function of switchings mean 
rate. In the final section we draw the conclusions. 



II. GENERAL EQUATIONS 

We consider the one-dimensional overdamped Brown- 
ian motion in switching potential profile 



dx 



d<^ (x,t) 
dx 



at). 



$ {x, t) = U (x) + V{x)r] (t) . 



(1) 



Here x (t) is the Brownian particle displacement, ^(t) is 
the white Gaussian noise with zero mean and correlation 
function (t) ( {t + t)) = 2DS (r). The potential $ {x, t) 
is the sum of two terms: the fixed potential U {x) and 
the randomly switching term V {x) rj (t) . The variable 
r] (t) is the Markovian dichotomous noise, which takes 
the values ±1 with the mean flipping rate v. If we invoke 
the following expression for probability density in terms 
of the average 



W{x,t) = {S {x-x{t))) 
and introduce auxiliary function Q {x,t) 

Q{x,t)^{f1 it)Six-x{t))), 



(2) 



(3) 



we obtain the next closed set of equations (see [23, |2€ 
and Appendix) 



dW d 

— ^-[U'ix)W + rix)Q] 



D- 



do 8 

-^ = -2.Q + -[U^ix)Q + V'ix)W]+D^^^ 



dx'^ 



(4) 



Let X — xq he the initial position of Brownian particles. 
Then 



W{x,0) ^Six-xo), 



(5) 



and W {x,t) becomes the conditional probability density 
W {x,t\xQ,0). Since ?y (0) is a deterministic value, the 
initial condition for the function Q {x, t) is (see Eqs. Q 
and ©) 



' (x, Q) = W [x, 0) ?7 (0) ^±5{x- xo) ■ 



(6) 



Let us consider the potential profiles U {x) ± V {x) with 
a wall at X ^ — oo and a sink at a: ^ -|-oo (see Fig. 
The potential profile U {x) + V {x) corresponds to a 
metastable state, and U [x) — V [x) corresponds to an 
unstable one. Thus, we investigate the system with ran- 
domly switching metastable state. 

The nonlinear relaxation time (NLRT) for the state 
located in the interval (Li, L2) is defined as follows |0| 



T (a^o) 



dt 



W {x,t\xo,Q)dx, (7) 



\U{x)-V{x) 




U{x)+V{x) 



FIG. 1; Switching potential with metastable state. 



where xq G (Li,L2). The NLRT is also interpreted 
as mean lifetime of Brownian particles in the interval 
(ii,L2) or average residence time, because, in accor- 
dance with Eq. 0, Eq. (fT)) can be rewritten as condi- 
tional time average 



' (x (t) - Li) e{L2-x [t)) dt \x (0) = xa 



where 9 [x) is the step function. 

Let us rewrite the definition Q in the form 



(xo) = / F(x,xo,0)dx, 



(8) 



where Y (x, xq, s) is the Laplace transform of conditional 
probability density VF (x, t |xo, 0). After Laplace trans- 
forming Eqs. I@J, with initial conditions and we 
obtain the following closed set of ordinary differential 
equations 

DY" + [U' (x) Y + V' (x) i?]' - sY = -5{x- xq) , 
DR" + [U'R + V'Y\ - (s + 2v)R = ip^ (x - xq) , (9) 

where i?(x,xo,s) is the Laplace transform of auxil- 
iary function Q{x,t), defined by Eq. ®. Using the 
method proposed in Ref. j27j |. we expand the functions 
sY (x, Xq, s) and sR (x, xq, s) in power series in s 

sY (x, xo, s) — Zq (x, xo) + sZi (x, xq) + . . . 

sR (x, Xo, s) = Rq (x, Xo) -I- sRi (x, xo) -I- . . . (10) 

Since all Brownian particles move to the sink located at 
the point x = +00 (see Fig. ^ we have zero stationary 
probability distribution, i.e. 

lim W (x, t |xo, 0) = lim sY (x, xo, s) = 0. 

t — *oo s — ^0 

As a consequence, in expansions (IIUII Zo(x,xo) — 0, 
Rq (x,xo) ~ 0, and the definition ^ becomes 



t(xo) 



Zi (x, Xo) dx. 



(11) 



Substituting the expansions (|10|l in Eqs. ® and equat- 
ing the terms without s, we obtain the following set of 
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equations for the functions Z\ (x, xq) and i?i (a;, a;o) 

DZ'l + [C/' (x) Zi + V^' (x) = -(5 (x - xo) , 

L'i?'/ + \U'Rx + V^'^i]' - ^vRx =^6(x~ xo) . (12) 

Because of the reflecting boundary at x = — cx3, the 
probability flow equals zero at this point, and from 
Eqs. Q we have 



dW 

D— + U' (x) W + V (x) Q 
+ U' (x) Q + V (x) 

ax 



= 0, 
= 0. 



(13) 



Making the Laplace transform of Eqs. p3|) and substi- 
tuting the expansions H1U|) . we obtain the following con- 
ditions for the functions Zi (x, xq) and Ri (x, xq) 



[DZ[ + U'ix)Z,+V' (a:)i?i],_^ = 0, 
[DR[ + U' (x) i?i + V (x) = 0. 



(14) 



By integrating the system (|12|l from — oo to x, with 
boundary conditions (|14|l , we obtain the following closed 
set of integro-differential equations for the functions 
Zi (x, xo) and Ri (x, xq) 

DZ[ + U' (x) Zi + V' (x) Ri = -6* (x - xo) , 
DR[ + U'Ri + V'Zi = 2v I Ridy T 6* (x - Xq) . (15) 

These general equations allow to calculate the NLRT for 
potential profiles above-defined. We may consider two 
mean lifetimes t+ (xq) and t_ (xq), depending on the 
initial configuration of the randomly switching potential 
profile $(x,0): U{x) + V{x) orU{x)-V{x). The NLRT 
(|ll|l is equal to r+ (xq), when we take the sign "— " in 
the second equation of system (|15|) . and vice versa for 

T_ (xq). 



III. LIFETIMES FOR PIECE- WISE LINEAR 
POTENTIAL 

Let us consider a piece- wise linear potential profile (see 
Fig. |2Jl with V {x) = ax {x > 0, < a < k) and 



-l-oo, X < 

U{x) = <( 0, < X < L 

k {L — x) , X > L 



(16) 



Hereafter we shall analyze the average residence time 
r_ (0) from the interval (Li = 0, L2 = b) with b > L, 
which is finite in deterministic case. We consider the 
initial position of all Brownian particles at the origin, i.e. 
Xq = 0. The potential profile U (x) -I- ax corresponds to 
metastable state and U (x) — ax corresponds to unstable 



A 


. U{x)+ax 
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FIG. 2: Switching piece-wise linear potential. 

one, as indicated in Fig. El After substituting the poten- 
tial H16|l and V (x) = ax in Eqs. H15|) and choosing the 
sign we arrive at 

DZ{ - k ■ 9{x - L) Zi+aRi = -1, 

DR[ - k ■ 9{x - L) Ri + aZi = l + 2iy [ Ridy. (17) 

We solve the set of differential equations H17|) in the re- 
gions < X < L and x > L separately, and then use the 
continuity conditions at the point x = L 



Zi |l-o — Zi |l+o , ^1 |l-o — Ri \l+o ■ 
For < X < L, the solutions of Eqs. (|17|l read 



(18) 



Zi (x) ~ ci [ cosh7x H — I -I- C2 smh7x H — 



a -f^D^ ' 



Ri (x) 



(ci sinh 7X -|- C2 cosh 7x) — — ■ . 

a 7^/9^ 



(19) 



where Zi (x) = Zi (x,0), i?i (x) = i?i (x,0) and 



a2 2v 



7 = \/ — + -■ (20) 
The finite solutions of Eqs. H17|l in the interval {L, +00) 



are 



Z, (x) =C3-e^(^-^) + |, 
k 



(21) 



where 



2k 
3D 



cos f = — 



^2£)2 

1 + 9 {lyD- a2) /p 



27r 



[l + 372i:)2/fc2]3/2 

is the negative root of the following cubic equation 



, , , k\ o , 2i'k 



(22) 



(23) 
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Substituting the solutions H19() and (|21|l in the continuity 
conditions ((TH)l and in the second equation ((T7|) . we ob- 
tain, after rearrangements, the following compact system 
of algebraic equations for unknown constants ci, C2, C3 

f2vk \ 

ci cosh7L + C2 smh7L + C3 — ^1 =0: 



/ir2 

• 1 r IT- k — llD 

d smh 7L + C2 cosh 7L + C3 — - — = - , (24) 



ci - C3 



where 



r = 7^, 

a 2vaL 
k 7^ 



/i = 



- 1. 



(25) 



The solutions of Eqs. are 



Cl 



C2 



C3 



a/i a^fc {2vDa sinh 7L - /iF^) 

2^ ^ 2iyr3i:) [a^fc + D {2vk - iiT^) cosh-/ L+iiTD [iiD - k) sinh7L] ' 
a{(^r2 - 2vk) {2vDa + hV^smh-^L) + [h^iV^ (k - ^jlD) - 2va^k] cosh7L} 
2i^r3 [a2fc + D {2vk - ^iV^) cosh'yL+fiTD {fiD - k) sinh7L] 
(2iyDasmh'-fL - hT^) 
2vT [a'^k + D {2vk - /iF^) cosh7L+/xrL> {jiD - k) sinh7L] ' 



(26) 



r 



Substituting the expressions 1 19|l and H21|) of the function 
Zi (x) in Eq. Hll() . and using Eqs. (|26|l , we obtain finally 



the following result for mean lifetime 



, , b uL^ 



DT {2va (e^'ib-L) _ 1) + 2vkL\ + hV'^ {2vk - ^F^)} sinh7L - 2vDa {a^k + ^iV'^D - 2vkD) 



2vT'^ ^ a2fc + D(2j/A:-^r2) cosh 7L+A*rD (/^D - fc) sinh 7L 

D [Vr* {ijlD -k)+ 2va {a^k + iiV'^D - 2vkD)\ cosh-fL - hV^ [r^ (e^(''--^) - l) + 2vkL + ijlD{hD - fc)] 

a2fc + D {2uk - /ir2) cosh7L+/irL» (/.tZ? - /c) sinh7L 

I 



)• (27) 



Equation H27|l is exact, and was derived without any obtained from Eqs. ()20|l and (|22|l . we find from Eq. H27|l 
assumptions on the white noise intensity D and on the 



mean rate of flippings v. 



r_(0) 



2 r 



h L 
k^2D 



1 



bq (1 - q) 
ujL 



(28) 



IV. CONDITION TO OBSERVE NOISE 
ENHANCED STABILITY (NES) 

Because of the complicated expression (|27|l of the mean 
lifetime, we analyze the limiting cases of very large and 
very small noise intensities. Using the approximate esti- 
mations for small parameters 7 and fi in the limit D —^ 00 



D \ AvD 



'-(1- ^ 



Here uj = vL/k and q = a/k are dimensionless parame- 
ters. The parameter q quantifies the degree of potential 
flatness after the point L (see Fig. I^J. Under very large 
noise intensity D, Brownian particles "do not see" the 
fine structure of potential profile and move as in the fixed 
potential —kx. Therefore NLRT decreases with noise in- 
tensity, tending to the value h/k as follows from Eq. H28|l. 

The NES phenomenon should be searched in opposite 
limiting case of very slow diffusion {D ^ Q) dill HI. 
The approximate expressions, obtained in this limit, for 
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parameters 7,r and /i are 
I'D 



7 



D 



1 



2uj 



r ~ a 1 



L(l-g2) 



2ujD (l + g2 
kL{l-q^f 



(29) 



terms up to first order in D, we obtain the following 
expression for NLRT at small noise intensity 



T_ (0) = TO + ^ • /(g, ^, s) + o . (30) 



Substituting Eqs. I|29|l in Eq. I|27|l . and retaining the Here 



fii^^^^)^—, — — + ^^— — — + j^-^^ — — + {^-^ ) — — — — (31) 



2(1-92) 



(l + g)(l-g2) 



2(1 -g2) 



and 



a a- 



is the mean lifetime in the absence of white Gaussian 
noise (D = 0). In Eqs. (|31|l and (|32|) new dimensionless 
parameter 

_ 2uj{h/L-l) 

is introduced. 

For very slow switchings — > 0, we find from Eq. H32|) 



2L b-L 
ro(i/^0) = — + -—, 
a k + a 

which is different from the deterministic time 

Td = To(iy = 0) = — + — — . 

a k + a 



(33) 



(34) 



Difference between the results (|33|) and (|34ll is due to 
a nonzero probability of one switching within the deter- 
ministic time interval {0,Td) in the case ^' — > 0. 

The condition to observe the NES effect can be ex- 
pressed by the inequality 



fiq,uj,s) > 0. 



(35) 



Let us analyze the structure of NES region on the plane 
(g, oj) from Eqs. H31|) and H35|) . At very slow and fast 
flippings we obtain 



9 > ^^^^^ 0,7863, 



to < 



qjSq^ + q-S) 
1-^2 



(36) 



In Fig. 13 we show the NES region (shaded area) on the 
plane {q,w) for b/L — 2, obtained from inequality (|35|l . 




FIG. 3: Shaded area is the parameter region on the plane 
(g,aj) where NES effect can be observed. Here u — {vL)/k, 
q = Q./fc, and b — 2L. 



The NES effect occurs at g ~ 1, i.e. at very small steep- 
ness k — a = k {1 — q) oi the reverse potential barrier for 
the metastable state. For this potential profile, a small 
noise intensity can return particles into potential well, 
after they crossed the point L. Then Brownian particles 
stay for long time in the metastable state. This means 
that, for a fixed mean flipping rate, the NES effect in- 
creases when g — !■ 1. For fixed parameter q the effect 
increases when oj ^ 0, because Brownian particles have 
enough time to move back into potential well. 

In Fig. ^ we show the plots of the normalized mean 
lifetime t_ (0) /tq, Eq. (|77|l . as a function of the noise 
intensity D for three values of the dimensionless mean 
flipping rate uj = vL/k: 0.01, 0.05, 0.1. The maximum 
value of the NLRT and the range of noise intensity values, 
where NES effect occurs, increases when uj decreases. 

By using exact Eq. (|27|l we have also investigated the 
behavior of the mean lifetime r_ (0) as a function of 
switchings mean rate v. In Fig. we plot this behavior 
for seven values of noise intensity. At very slow flippings 
[u —t 0) we obtain 



(0) ^ Td 



D \ 



-aL/D\ 



a2 (1 + q) 



(37) 
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FIG. 4: Semilogarithmic plot of the normalized mean lifetime 
r_ (0) /ro vs the white noise intensity D for three values of 
the dimensionless mean flipping rate lj — vL/k: 0.1 (curve 1), 
0.05 (curve 2), 0.01 (curve 3). Parameters are L = 1, fc = 1, 
6 = 2, and o = 0.995. 
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This suggests that, the enhancement of stability of 
metastable state is strongly correlated with the poten- 
tial fluctuations, when the Brownian particle "sees" the 
barrier of the metastable state 



)wnian par t 



V. CONCLUSIONS 

We have investigated the nonlinear relaxation time 
for one-dimensional system with additive white Gaussian 
noise, and potential profile switching between two config- 
urations, due to a Markovian dichotomous noise. From 
the general equations H15|l . we provide exact expression 
of the mean lifetime for piece- wise linear potential, for ar- 
bitrary noise intensity, and arbitrary fluctuation rate of 
the potential. We find the noise enhanced stability and 
the resonant activation phenomena in the system investi- 
gated. We obtained analytically the region on the (g, lu) 
plane, where the NES effect can be observed. Moreover, 
when we fix white noise intensity D, flatness q, and vary 
switchings mean rate v, we can observe new resonant- 
like behavior of the mean lifetime, which is related to the 
NES phenomenon. The NLRT shows a maximum as a 
function of the mean flipping rate of potential, with the 
NES effect strongly correlated with the potential fluctua- 
tions. The general equations derived in this paper enable 
us to perform the analysis of the NES effect conditions in 
physical systems with more complex potential profiles. 



FIG. 5: Semilogarithmic plot of the mean lifetime r_ (0) vs 
the dimensionless mean flipping rate u = vL / k for seven noise 
intensity values. Specifically from top to bottom on the right 
side of the figure: D = 0.03,0.05,0.07,0.09,0.15,0.25,0.35. 
The other parameters are the same as in Fig. 0| 



i.e. the NLRT of the fixed unstable potential U [x) — ax. 
While for very fast switchings [v — + oo) we obtain 



r_(0) 



2D' 



(38) 



i.e. the mean lifetime for average potential U [x). All 
limiting values of the NLRT expressed by Eqs. (|37|l and 
(|38|l are shown in Fig. El At intermediate rates the escape 
from the metastable state exhibits a minimum at w = 
0.1, which is the signature of resonant activation (RA) 
phenomenon [3. H. Ha. . 

Moreover, in Fig. [S] we observe a new resonant-like 
behavior for the NLRT as a function of mean fluctua- 
tion rate of potential. The NLRT exhibits a maximum 
between the slow limit of potential fluctuations (static 
limit) and the RA minimum. This maximum occurs for 
a value of the barrier fluctuation rate on the order of the 
inverse of the time T„p [D] required to escape from the 
metastable fixed configuration 



b-L L D (e°-^/^ - 1) 



a2 (1 - q) 



(39) 
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APPENDIX 

Upon differentiation of Eq. (01 on t, we obtain 
dW d 

-^^-^{i{t)5{x-xm- (40) 

Substituting x it) from Eq. , and using the definition 
(|^ of auxiliary function, we can rewrite Eq. (|40|l as 



dW 
~dt 



d_ 

dx 
d_ 

dx 



[U'{x)W] 



d_ 

dx 



[V'{x)Q] 



{at)Six-x{m 



(41) 



To obtain the evolution of function Q(x,t), we use the 
Shapiro-Loginov's formula for Markovian dichotomous 
noise ^3 



dt 



{vit)Rt[rj]) ^-2iy(r^{t)Rt[fi]) + {v{t)Rt[v]), (42) 
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where Rt [ij] is an arbitrary functional depending on the 
history of random process f] (t), < t < t. Replacing 
Rt[r]] with S{x-x{t)) in Eq. (g2l, using Eq. (QJ and 
taking into account that ij^ (i) = 1, we arrive at 



dQ 
dt 



~2vQ 



d_ 

dx 



[U'{x)Q] 



d_ 

dx 



[V'{x)W] 



dx 



{at)ri{t)5{x~x{t))). 



(43) 



To split the functional averages in Eqs. H4HI and H43|l . 
we use the Furutsu-Novikov's formula for white Gaussian 
noise ^ {t) JflJ 



dT 



D 



sat) 



(44) 



where Ft[£] is an arbitrary functional of £,{t). Re- 
placing sequentially i^t [C] with S {x — x (t)) and with 
7] (t) S {x — X {t)) in Eq. H44|l . and taking into account 
that, in accordance with Eq. (Q, Sx {t) /5£,{t) — 1, we 
find 



{^{t)6{x-x{t))) = -D 
{£,{t)r^{t)5{x~xm^-D 



dW 
dx 
dQ 
dx ' 



(45) 



Substituting the expressions (|45|l in Eqs. (|41|1 and H43|l. 
we obtain the desired closed set of equations 
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